Spectrum of coherently backscattered light from two atoms 
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We present a detailed analytical and numerical analysis of the inelastic coherent backscattering 
spectrum of laser light incident on cold atoms. We identify frequency domains where the interference 
contribution can be positive as well as negative - or exhibits dispersive character. These distinctive 
features are explained by reciprocity arguments and dressed state two-photon scattering amplitudes. 

PACS numbers: 



I. INTRODUCTION 

Multiple scattering of light in cold atomic gases has be- 
come an area of intense theoretical and experimental re- 
search (for a recent review, see [l[). On the experimental 
side, the successful observation of coherent backscatter- 
ing (CBS) of light in clouds of cold atoms 0, Q demon- 
strated the potential of finely tunable atomic media for 
detailed studies of localization and transport phenomena 
4] in the weak and, prospectively, strong localization 
regime. On the theoretical side, it is of crucial impor- 
tance to understand how interference effects are affected 
by various dephasing mechanisms characteristic for atom- 
photon interactions, for ensembles of atoms cooled down 
to approximately 100 /xK, with, in general, degenerate 
electronic structure, and in the presence of inelastic scat- 
tering. It is presently also realized that multiple scatter- 
ing of light in atomic clouds is relevant in the context of 
quantum information storage and retrieval by photons, 
in the parameter regime of electromagnetically-induced 
transparency (EIT) 0. This promotes CBS studies to- 
ward the dynamical regime [6(. Also multiple scattering 
of nonclassical light Q is presently moving into focus. 

In our present contribution, we will expand on the im- 
pact of inelastic scattering processes on CBS. Experimen- 
tal studies on cold Sr atoms revealed a rapid decrease of 
the CBS interference contrast with increasing intensity of 
the injected laser field, as a consequence of the saturation 
of the laser-driven atomic transition [8| . A much weaker 
sensitivity of the quality of the CBS signal was observed 
for Rb atoms [9(. In this latter case, inelastic processes 
occur on degenerate transitions [HI EH, E3, EE E3] > while 
Sr atoms offer dipole transitions with a nondegenerate 
ground state [l5| . 

So far, the role of the nonlinearity of the atom-photon 
interaction for CBS scatterers has been theoretically in- 
vestigated only for Sr atoms 0, OJ, El EL E| . Within 
a scattering theoretical approach applied to two atoms 
in the regime of weakly nonlinear scattering (l6j . the 
decrease of the CBS enhancement factor a was shown 
to be due to the partial distinguishability of the inter- 
fering amplitudes. In the general case of many atoms, 



three different amplitudes interfere constructively in the 
weakly nonlinear re gim e, such that a may exceed the lin- 
ear barrier two fl7l Il8|| . For arbitrary intensities of the 
injected laser light, a number of effects have been pre- 
dicted within a master equation approach |19|,|20( , such as 
a nonvanishing residual CBS contrast in the deep satura- 
tion regime, or CBS anti-enhancement under off-resonant 
driving. The strongly inelastic scattering regime studied 
by a quantum Langevin treatment highlighted the crucial 
role of inelastic susceptibilities [2lj . 

Here, we will focus on the spectral properties of the 
CBS signal, in order to provide a detailed interpretation 
of the residual CBS enhancement or anti-enhancement 
predicted earlier 0, El- This will also elucidate the 
structure of the CBS spectra presented in ELEI ■ Specif- 
ically, we will identify frequency domains where the inter- 
ference contribution to CBS exhibits not only construc- 
tive but also destructive terms, or else a dispersive line- 
shape. We relate the spectral lines of the CBS spectrum 
to CBS transitions between atomic dressed states. We 
show that, in the limit of very intense laser fields, spec- 
tral ranges with destructive interference lead to a de- 
creased enhancement factor, which, at exact resonance, 
shrinks to — 23/21. For off- resonance driving, the 
destructively interfering processes can outweigh the con- 
structively interfering ones, leading to values of the en- 
hancement factor less than unity. 

The paper is organized as follows: We start with a brief 
presentation of our two-atom model, and of the master 
equation approach. In Sec. IIIII we present our results on 
the CBS spectrum. Section [TVl concludes the paper. 



II. A MASTER EQUATION FOR TWO ATOMS 

A. The model, and the main quantity of interest 

While details of our approach were presented elsewhere 
[20l ]. we recollect its basic ingredients relevant for our 
subsequent spectral analysis. We start out with the gen- 
eral formulation of a Hamiltonian describing N identical, 
stationary atoms embedded in an electromagnetic envi- 
ronment of quantized harmonic oscillators, and subjected 
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FIG. 1: Model of CBS with two atoms, (a) atoms (black 
dots) are driven by laser light with right circular polariza- 
tion, while CBS is observed in the helicity preserving channel, 
that is, with flipped polarization. Photons in this channel ap- 
pear as a result of double scattering, g is the strength of the 
far-field dipole-dipole coupling responsible for the exchange 
of photons; (b) internal atomic structure corresponding to a 
J g = — > J e = 1 dipole transition. u>o is the transition fre- 
quency, 27 is the radiative linewidth, Q is the Rabi frequency. 
Sublevels |1) and |3) have magnetic quantum number m = 0. 
Sublevels \2) and |4) correspond to m = — 1 and m = 1, re- 
spectively. The thick solid arrow represents the laser field 
driving the |1) <-> |4) transition, while the dashed arrow indi- 
cates the CBS field originating from the |1) <-> |2) transition. 



to an external (classical) laser field of arbitrary intensity. 
Coupling to the bath gives rise to spontaneous emission 
from the excited state and to the far-field dipole-dipole 
interaction responsible for the exchange of photons be- 
tween the atoms, whereas the coupling to the laser field 
induces Rabi oscillations of the populations and coher- 
ences on the laser-driven atomic transitions. 

The intensity of the light scattered off the atomic sys- 
tem is expressed via the correlation functions of the emit- 
ting dipoles. We address a regime of multiple scattering 
which is relatively simple from the theoretical point of 
view: scattering by an optically thin atomic medium, 
where double scattering provides the dominant contri- 
bution to the CBS signal [1(1 EH- It is in this double 
scattering regime where the first observation of a CBS 
reduction due to the saturation of atomic dipole transi- 
tions was reported Q. 

Under this specific conditions, considering a system of 
only two atoms suffices to grasp the essential physical 
phenomena. Thus, we come up with the toy model of 
CBS depicted in Fig. [Q We will study CBS from two 
identical, motionless atoms, located at positions rj and 
r2, with the distance r\i = \r± — r%\ being much greater 
than the optical wavelength A. 

As for the internal atomic structure we choose non- 
degenerate atomic ground states, and the excited state 
with a three- fold degeneracy [see Fig. Hfjb)], precisely as 
in the Sr experiment Q. The laser intensity is encoded 
in the saturation parameter s = £l 2 /2(~f 2 + 6 2 ), where 

is the Rabi frequency, 7 is half the spontaneous decay 
rate of the atomic excited states, and S = lol — ^0 is the 
laser-atom detuning. 

We will consider the CBS signal in the helicity pre- 
serving (h || h) polarization channel, as in [8j|, with right 



circularly polarized laser light driving the |1) «-> |4) tran- 
sition, that is, £l = e + i, in helicity basis notation. CBS 
with preserved helicity then corresponds to the detection 
of photons with flipped polarization, e = e_i, that is, 
from |1) <-* |2) transition, as shown in Fig.[TJ 

The CBS spectrum can be derived from the average 
value of the first-order temporal correlation function of 
the field 

G« (r, t; r,t') = /Tt {p[e ■ E<"> (r, t)] [e* ■ E« (r, f )]}) 

\ / cor 

where p is the initial density operator of the atom-field 
system, Es (r,i) is the negative/positive frequency 
component of the electric field operator of the scattered 
field, and (. . .) conf denotes a configuration average. The 
components of the scattered field are the retarded fields 
radiated by the atomic dipoles, 

E(+)(r,<) = ^i^^D a (t a )e- <kr <» , (2) 

Q — 1 

where e is the permittivity of the vacuum, D a = 
-e-iof 2 + e of 3 - e + icrf 4 , with a% l = \k) a (l\ a , is the 
dipole lowering operator, and t a = £ — |r — r a \/c. In 
writing Eq. |2|). we have assumed that r\i <ti r, that 
is, the field is detected at a distance much larger than 
the interatomic distance. In the following, we will for 
brevity omit the r-dependent prefactor of Eq. @ and, 
consistently, of the temporal correlation functions. 

Inserting Eq. © into Eq. (|TJ) we obtain, in the steady 
state limit t — ► 00, 

2 

G${t) = ((^i^ 2 (r)) ss e ik - r ^)conf, (3) 

where "ss" stands for steady state, r = t' — t > 0, 
r Q/ 3 = r Q — rp, and the inner angular brackets indicate 
the quantum mechanical expectation value [see Eq. (fTj)]. 
The spectrum follows via a Laplace transform of ^ 

^) = ~lim Re{G«(z)}, (4) 

where Gas (z) = J™ drexp(-zr)Gss ) (r), z = T — iu, 
with r > 0, and v = lo — u>l. Note that the spectrum is 
defined with respect to the laser frequency, what implies 
that the atomic correlation functions must be evaluated 
in the frame rotating at ljl- 

B. Configuration average 

A configuration average is necessary because the two- 
atom correlation functions may sensitively depend on the 
interatomic distance, and on the orientation of the vector 
ri2 with respect to k^, and thus exhibit rapid oscillations 
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around the backscattering direction. These oscillations 
have the same nature as a speckle pattern scattered off a 
disordered medium. After many realizations of the dis- 
order all interference maxima except the one due to CBS 
disappear. A simple and sufficient way to mimic disorder 
in our two-atom system is to assume an isotropic distri- 
bution of the radius-vector connecting the atoms, and a 
uniform distribution (with a width A) of interatomic dis- 
tances around the average distance £ equal to the scat- 
tering mean free path. 



C. Master equation 

To deduce the atomic correlation functions which en- 
ter the right hand side of Eg. ((3J we adapted [H, H3| 
the theoretical approach of [25| . Within this setting, the 
dynamics of the dipole operators' expectation values as 
well as of the dipole-dipole correlators is governed by the 
master equation 



(5) 



where the Liouvillians C a and C a p generate the time evo- 
lution of an arbitrary atomic operator Q, for independent 
and interacting atoms, respectively. Explicitly, 

C a Q = -i5[T>t ■ D a , Q] - l - [n Q (D+ ■ e L ) + Q* a (D a • e* L ) , Q] 

+7(D^[Q,Dj + [Dt,Q].D a ), (6) 
C aP Q =Dt -T(g, n) • [Q, Dp] + [D+ , Q] ■ T*(g, n) D Q ,(7) 

where f2 Q = f2e lkl, ' ra is the position-dependent Rabi fre- 
quency. The radiative dipole-dipole interaction due to 
exchange of photons between the atoms is described by 
the tensor T (g, n) = -fg A , with A = 1 — nn the pro- 
jector on the transverse plane defined by the unit vector 
n along the connecting line between atoms a and (3. This 
interaction has a certain strength, depending on the dis- 
tance between the atoms, via 
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2k r af: 



(8) 



with fco = wo/c, and on the life time of the excited atomic 
levels, through 7. The coupling constant \g\ <C 1 is small 
in the far-field (fco?"a/3 3> 1), where near- field interaction 
terms of order {k^r a p)~ 2 and (kor a p)~ 3 can be neglected. 

Of course, an arbitrary operator Q inserted into Eq. ([5]) 
does not result in a closed differential equation. Our 
system consisting of two 4-level atoms leads to 255 = 
4 2 x 4 2 — 1 linear coupled equations of motion for the 
associated expectation values. We solve them pertur- 
batively to second order in g, to account for the lowest 
order (double-) scattering processes giving rise to a non- 
trivial interference contribution. To keep this in mind, 



symbols denoting double scattering intensities and spec- 
tra will carry the subscript 2. 

Note that Eq. describes the evolution of one-point 
correlation functions, whereas (r) is a two-point cor- 
relation function. By virtue of the quantum regression 
theorem (24|, the latter also satisfies Eq. ([5]), with ini- 
tial conditions extracted from the stationary solution of 
([5l). In particular, the double scattering counterpart of 

Gil\o) is nothing but the stationary average backscat- 
tered light intensity which will be referred to as I^- 
There is an obvious relation between and 5a (f) [ob- 
tained by expanding © to second order, in Q]: 



rtot 

1 2 — 



(9) 



The total CBS intensity at backscattering direction can 
be decomposed in a sum of two terms, 



jtot 7^°* 



Ctot 
2 ! 



0) (i.e., k 



and 



LT = ((^2>I 2 S 1 + (^2> [2l >ss>co„f, 



(10) 

(11) 
(12) 



are the so-called "crossed" {C\ ot (6 — 0)) and "ladder" 
(£2°*) terms, respectively. Using these, we can derive the 
standard measure of the phase coherence between the 
counterpropagating amplitudes that contribute to CBS 
- the enhancement factor 



1 



T tot " 

^2 



(13) 



For perfect two- wave interference, a — 2. In general, the 
total backscattered light intensity has elastic and inelas- 
tic components, 



rtot rel 

1 2 — L 2 



(14) 



with the elastic component given by products of the ex- 
pectation values of the atomic dipolcs, 



4= E ((^{^e 



ikr,-r Q 



(15) 



a, 0=1 



For a = (3 we obtain the elastic ladder term Lf, and for 
a (3 the elastic crossed term Cf . Given I| ot and if , 
also the fluctuating part of the dipole correlation func- 
tions defining 7 2 nel is determined. 

A detailed derivation of the stationary CBS intensity, 
together with analytical results for S — 0, can be found 
in Appendix A. 



III. CBS SPECTRUM 

With the above premises, we can now proceed to the 
detailed analysis of the CBS spectrum, the central object 
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FIG. 2: (Color online) Elastic ladder, L| (solid), inelastic 
ladder, L 1 ^ 1 (dashed), and inelastic crossed, Cjj , (dashed- 
dotted) CBS contributions vs. Rabi frequency Q, for different 
values of the detuning 8: 8 = (a) 0; (b) 57; (c) 2O7. 



of this paper. A detailed derivation of the CBS spectrum 
from the solution of Eq. ([5]) is given in the Appendix B, 
while the physical content thereof will be discussed in the 
following sections. 



A. Elastic spectrum 

In the saturation regime, the double scattering spec- 
trum of CBS has elastic and inelastic parts. The elastic 
spectrum in the backscattering direction reads 



(16) 



where 8{v) is Dirac's delta- function, and if = Lf + Cf , 
with the ladder and crossed contributions [2fJ 

2\g\ 2 1 



TCI f~1Cl 

Ll 2 — °2 — 



15 i + (<y 7 ) 2 (i + s y 



(17) 



where g is the configuration averaged value of g. A de- 
tailed interpretation of the elastic intensity was given in 
previous work [l9l [20|. Here, we will focus on the . . . 
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FIG. 3: (Color online) Normalized inelastic spectra of the 
ladder, U^iv) (solid), and crossed, C^^lv) (dashed), terms 
at n = O.I7. 8 = (a) 0; (b) 7; (c) 5 7 . 



the asymptotic value of the interference contrast of CBS, 
Ooo(^) — 1- Furthermore, for arbitrary saturation, the 
areas under the peaks of the normalized ladder spectrum 
give the relative probabilities of the corresponding inelas- 
tic processes. 

The normalization itself depends on the parameters of 
the driving field. In Fig. [2j we present several examples of 
the elastic and inelastic intensities, as functions of 17, and 
for different values of the detuning S. This figure shows, 
in particular, that, for <5 = 2O7, the inelastic crossed 
term is negative, for a range of Rabi frequencies, with a 
minimum at £1 ~ 2O7 = 5 (see Fig. [H[c)). Furthermore, 
note that, since Cf 1 = L| , the total crossed term Cf 1 
1 < around Q, = 2O7 
and 4 neI 

2.82 x 10~ 5 , and 4.27 x 10 -4 , respectively, what implies 



The values of Cf = ~L$, 
L^ cl at n = 2O7 and S = 2O7 are 2.45 x 1CT 5 , 



C^ el /Lf ei = -0.066, and a = 0.991. The negativity of 
the total crossed term thus results in anti-enhancement 
(i.e., an enhancement factor a < 1) and was reported 
previously [20]. It is one of the purposes of the present 
paper to identify the physical origin thereof. 



2. Weakly inelastic scattering 



B. Inelastic spectrum 

1. Normalization 

Since we are interested in the spectral properties of 
CBS, we will omit all frequency-independent prefactors. 
This is naturally achieved when considering normalized 
expressions. We choose the stationary inelastic ladder 
contribution L™ el as normalization factor, such that the 
integrals of Lf c \v)/Lf cl and C^ c \v)/Lf cl over v yield 
unity and C^/L™ 1 , respectively. In the deep satu- 
ration regime, the value of the latter integral tends to 



In the regime of small Rabi frequencies fl -C 7, it is the 
lowest-order inelastic processes - two-photon processes, 
proportional to O 2 - that contribute to the inelastic spec- 
trum. Consequently, the emerging spectral features of 
the ladder and crossed spectra derived with the aid of 
the master equation approach (see Fig. [3]) can be in- 
terpreted on the basis of two-photon amplitudes derived 
in [161 ] . There it was shown that, in the weakly inelas- 
tic regime, one of the two atoms must scatter inelasti- 
cally. Furthermore, for direct and (time-)reversed am- 
plitudes (which have to interfere constructively to create 
the CBS signal), the inelastically scattering atom must 
be the same. Then it turns out that, although initial and 
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(a) 
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FIG. 4: Direct (a) and (time-)reversed (b) scattering 
paths/amplitudes, giving rise to the CBS signal in the weakly 
inelastic scattering regime. Grey-shaded and black spots rep- 
resent inelastically and elastically scattering scattering atoms, 
respectively. The amplitudes spell out the transformation of 
ingoing frequencies oil into outgoing frequencies u>. These in- 
terfering amplitudes are non-reciprocal, since the frequencies 
of the intermediate photons for direct and reversed paths are 
different. 



final frequencies of the scattered photons are the same, 
the intermediate frequencies differ, leading to the non- 
reciprocity of the interfering amplitudes, and a decrease 
of the CBS enhancement factor. The non-reciprocity ar- 
gument was also successfully applied to interpret a dis- 
persive resonance of the crossed term's spectrum in the 
weakly inelastic regime [21] . 

Figure Q] shows direct (a) and reversed (d) processes 
contributing to the CBS signal. Without loss of gen- 
erality, we assume that the left atom is scattering in- 
elastically. The direct, E\, and reversed, E 2 , scattering 
amplitudes can be derived as [l6| 



Ei = 
E<x — 



(7- 



- i(8 - i/)) (7 

e i<t>/2 



i(8 + v)) s 



(7 + i5) 2 + v 2 



(18) 
(19) 



where is a phase which depends on the geometric con- 
figuration, and should be taken zero in the backscattering 
direction. From Eqs. (|18[l and (fl"9|) . we obtain following 
expressions for the ladder and crossed spectra: 



2( 7 2 + S 2 ) + 2Sv + v 2 
( 7 2 + (5 - is) 2 )(j 2 + {5 + u) 2 ) 2 
2{ 1 2 +5{5 + v)) 

( 7 2 + (5_ „)2)( 7 2 + ( ( J + I/ )2)2 



, (20) 
• (21) 



Expressions (|2"0j) and (|2ip precisely reproduce the spec- 
tral lineshape of the ladder and crossed spectra of Fig. [31 
and, up to the prefactor £1 4 /(87T7), coincide with the re- 
spective results of the master equation approach at van- 
ishing detuning 5 = (see Eqs. (|B13j) and (|B14j) in Ap- 
pendix B). They also allow for a transparent interpreta- 
tion of the CBS spectra for arbitrary detuning 5. 

By inspecting the denominators of Eqs. ([2"TJ)) and ([2~T) . 
we see that, at 5 = 0, both Lf e \v) and Cf el (v) must 
exhibit a single peak, at v = 0, as is the case in Fig.^a). 




4> |2> 



(a) § 




(b) 5 



FIG. 5: Two-photon processes describing the laser frequency 
transformation on the left atom in Fig. [4] (a) corresponds to 
Fig.[4ja): Two laser photons are converted into two scattered 
photons such that uj\ ~ ujq + 25 and u>2 — Uo- Either one 
of these photons then travels to the other atom, (b) corre- 
sponds to Fig. H|b) in case lu ~ ui : One laser photon and one 
elastically scattered photon give rise to the CBS photon at 
lu ~ u>i and the undetected fluorescence photon at LO2 — LUo. 
The roles of the photons at wi and u>2 can be interchanged by 
flipping their polarizations. 



At 5 =/= 0, there must be two resonances at v ~ ±<5, with 
the more pronounced one at v ~ —5 (since the respective 
term is squared) (see Fig. [3][b,c)). The physical reason 
for this behavior is easy to understand by recalling that 
upon inelastic scattering on the (left, in Fig. 2]) atom, the 
re-emitted photon frequency ui can be either ojl — 5 = ujq 
or u>l + 5 = u>q + 25, as spelled out by the level diagram 
in Fig. The former, equal to the atomic resonance 
frequency, corresponds to v ~ —5, and acquires a large 
scattering cross-section, while the latter corresponds to 
v ~ +5, and has a diminished scattering cross-section 
since detuned by 28 from the atomic transition frequency. 

Another peculiarity of the CBS spectra which clearly 
stands out in Fig. [3jc) as well as in Eq. ([21]) is a disper- 
sive lineshape of the crossed spectrum around v ~ —6. In 
other words, as the signal frequency passes from lu > ujq 
to w < ujq, the interference character between the in- 
terference paths (a) and (b) of Fig. 2] turns from con- 
structive to destructive. This is due to the continuous 
^-dependence of the phase shift between the direct (fTS]) 
and reversed (|19p scattering amplitudes, passing through 
7r/2 at lu ~ luq. 



3. Strongly inelastic scattering 

As the Rabi frequency fl increases to values ft > 7, 
scattering processes of higher than second order con- 
tribute for the individual atoms, resulting in the emission 
of the resonance fluorescence Mollow triplet [26j . Corre- 
spondingly, the CBS spectra become more complicated at 
intense driving. Figure [S] presents examples of normal- 
ized ladder and crossed spectra for two fixed values of 
the Rabi frequency, £1 = IO7 (Fig. [BJa-c)), and ft = 2O7 
(Fig.Htd-f)). 

At exact resonance, the spectra remain symmetric, 
with signatures of five resonances at v = ±f2, v = ±£7/2, 
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FIG. 6: (Color online) Normalized inelastic spectra of the 
ladder, L2 ncl , (solid) and crossed, C™ 1 , (dashed line) terms. 
Left panel ((a)-(c)): = IO7; right panel ((d)-(f)): Q = 20-y. 
The detuning increases from top to bottom: 5 = (top); 57 
(middle); 2O7 (bottom). 



and v — 0. Crossed spectra take negative values in a 
frequencies range which spares out the central, positive, 
peak. With increasing (5, some of the resonances ap- 
proach each other (observe, e. g., the change of position 
of the two left-most resonances of the ladder term from 
Fig. [fj(a) to (b), and from (d) to (e)), and eventually 
merge (see Fig- Et^f))- Fig- E(c) is qualitatively reminis- 
cent of Fig. [3](c) , though additionally garnished with a 
signal at v ~ in both the ladder and crossed spectra. 
This additional resonance stems from the central peak of 
the Mollow triplet, emerging in three-photon scattering 
processes on one atom. In contrast, Fig. [3fc) corresponds 
to a much weaker saturation parameter (s = 2 X 10~ 4 ) 
where three-photon processes are negligible, and hence 
the central Mollow peak is not visible.. 

A notable feature of the spectra of Fig. [6] is that, for 
5 — 2O7 ((c) and (f)), the overall inelastic crossed term 
becomes negative, i.e., f™^ duC^ e \iy)/L^ ei = -0.029 
(for Q, = IO7), and -0.066 (for (1 = 2O7), respec- 
tively. The negative value of the normalized crossed 
term in the latter case dominates the positive ratio 
Cf/Lf el = 0.057 which can be extracted from Fig. [3(c) 
for £1 = 2O7. Therefore, the enhancement factor a 
1 + (Cf + C$ iel )/(Lf + 4 nGl ) is e qual to 0.991 < 1, 
in accordance with [20j | and Sec. IHI B 11 It is now 
clear that the anti-enhancement comes from destructive 
self-interference of inelastically scattered photons around 

S 



v = -287 ~ -v 7 ^ 2 + S 2 , with n = 6 = 2O7. As seen 
in Fig. H(f), a t these parameter values, the two reso- 
nances at which the crossed spectrum is negative are 
almost merged. Hence, the destructively interfering pro- 
cesses responsible for each distinct resonance become in- 
distinguishable, and interfere with each other. This in- 
terference additionally broadens the negative area of the 
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FIG. 7: Single (left) and double (right) scattering processes. 
Only the laser-driven and the CBS transition are depicted, 
respectively. Radiative processes (a) and (b) on the left give 
rise to the Rayleigh peak of the Mollow triplet with linewidth 
7 centered at v = 0. (c) and (d) give rise to the sidebands 
at v = Q and v = —SI, respectively, with linewidths 37/2. 
Photons emitted by atom 1 (left) propagate to atom 2 (right), 
and are scattered from either one of the dressed states of the 
energy manifold with TV — 1 laser photons. The diagrams on 
the right describe the scattering of: (a)-(d) a sideband photon 
centered at v — — Q; (e)-(h) a photon from the central peak; 
(i)-(l) a sideband photon at v — Q; (m), (n) resonant photons 
at v = ±£1/2. 



crossed spectrum, leading to anti-enhancement. As the 
detuning becomes much larger than the Rabi frequency, 
only the dispersive resonance of the crossed spectrum sur- 
vives, leading to a picture qualitatively similar to that of 
Figs. OKc), [5(c). The position of this resonance corre- 
sponds to the resonance frequency of an ac-Stark shifted 
CBS transition, and will be specified in the next subsec- 
tion. 



4- Limit of well-separated spectral lines 

Although the positions of the resonances of the CBS 
spectrum could be guessed from the results presented in 
the previous section, the spectral line shape for the ladder 
and crossed terms is not quite clear, since not yet fully 
resolved at the Rabi frequency considered. To reach a 
fully transparent picture, we will now address the limit 
of well-separated spectral lines, at £1 ^> 7. 

a. Exact resonance. Let us start with a reminder of 
the level structure of the pumped (by the driving laser) 
and probed (by the scattered photon) transitions of both 
atoms. In an intense laser field with f2 3> 7, the pumped 
transition |1) «-> |4) is strongly coupled to the laser mode. 
In this case, it is instructive to treat the latter as a quan- 
tum system [13, [Hj]. The eigenstates of the laser-atom 
interaction Hamiltonian are the dressed states | ± (N)) a . 
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For 6 = 0, they read 

\±(N)) a = ^=(\l,N+l) a ±e ik - 



•|4,JV)a), (22) 



where iV and N + 1 refer to the number of photons in 
the laser mode, and a labels the atoms. Inasmuch as 
the dressed states represent superpositions of the ground 
and excited atomic states of the laser-driven transition, 
they have widths equal to 7 for the resonant driving [28| . 
Spontaneous transitions from the dressed states manifold 
{| ± {N)) a } to {| ± (N - 1)) Q } lead to emission of the 
fluorescence spectrum centered at frequencies wl- O, ui, 
and Wi+O known as the Mollow triplet [26| (see the left 
of Fig. [7]). The width of the central peak is defined by the 
decay rate of the dressed state's populations, and is equal 
to 27, while the width of the sidebands is determined 
by the decay rate of the coherences between the dressed 
states, and is equal to 37 [28| . 

On the other hand, the level |2) is not affected by the 
laser field. As a result, the probed transition |1) «-> |2) 
of both atoms, giving rise to CBS photons, is affected by 
the laser field only via the state |1). For that reason, new 
resonance frequencies of the probed transition emerge at 
uj l ± 0/2 (the right of Fig. 

When the Mollow triplet emitted by one atom is in- 
cident on another atom, it is scattered on the internal 
structure of the latter. Relevant scattering processes 
that can take place are depicted on the right of Fig. [7] 
Each photon can be scattered either elastically or un- 
dergo Raman-Stokes or -anti-Stokes (multiphoton) tran- 
sitions [29] (which lead to a frequency change by — O or 
O, respectively), which conserve energy and angular mo- 
mentum. It follows that the CBS spectrum must have 
resonances at v = ±20, v = ±0, and v = 0. The dia- 
gramms describing the emission of CBS photons at these 
frequencies are (a)-(l) on the right of Fig. [7] 

However, there are two more diagrams, (m) and (n), 
that apparently do not fit into the scheme just described. 
On these diagrams, photons with frequencies v = ±0/2 
are resonantly scattered to give rise to an additional dou- 
blet in the CBS spectrum. This is nothing but an Autler- 
Townes doublet [30 . |31| : One should bear in mind that 
a Lorentzian distribution has long tails, and thus emis- 
sion of photons with frequencies very different from the 
central frequencies is not impossible, though with very 
small probabilities. For instance, two laser photons may 
be transformed into two fluorescence photons with fre- 
quencies ujl + 0/2 and ujl — 0/2, with each of these 
photons appearing as a result of quantum interference 
between the Rayleigh and Raman-anti-Stokes or -Stokes 
transitions, respectively [Hj]. 

What makes the frequencies u>l ±0/2 special in our 
present problem is that both (dressed) atoms have two 
transitions exactly in resonance with these frequencies. 
Therefore, although the probability of the creation of a 
pair of photons with frequencies ± 0/2 is relatively 
small, the probability of their double scattering is high, 
due to a large resonant scattering cross-section on the 
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FIG. 8: (Color online) Normalized inelastic spectra of the 
ladder, Z' 2 ncl (i/) (solid), and crossed, C™ el (dashed), terms, in 
the limit of well separated spectral lines, at f2 = IOO7. 5 = (a) 
0; (b) 2O7. The numbers near the resonances indicate their 
areas, such that the overall areas of the ladder and crossed 
terms give unity and C™ 1 / L'^ el , respectively. 



transition |1) <-> |2) modified by the ac-Stark effect. 
The overall effect finally acquires the same strength as 
the elastic (though non-resonant) scattering of photons 
which are emitted into the dominant frequency compo- 
nents of the Mollow triplet, by the first atom. 

The results of our calculation are presented in 
Fig.[5Ja). Consistently with our analysis, both the ladder 
and crossed spectra consist of seven resonances, with the 
positions precisely at the predicted frequencies, that is, 
at v = ±20, v = ±0, v = ±0/2, v = 0. The respective 
analytic expressions, (|B19|) and (|B20|) . can be found in 
Appendix B. 

Let us now discuss the ladder spectrum in more detail, 
while postponing the analysis of the spectrum of the in- 
terference contribution for the next subsection. In the 
ladder spectrum, three of its resonances, at v = and 
v = ±0, represent sums of two Lorentzians with differ- 
ent widths. The remaining ones are simple Lorentzians. 
In order to understand this structure quantitatively, let 
us reinspect the diagrams on the right of Fig. [7] We 
start with the central CBS resonance at v = 0. Several 
processes, (d), (e), (g), and (j) on the right of Fig. [JJcon- 
tribute. Diagrams (e) and (g) represent a photon emit- 
ted by the first atom into the inelastic Rayleigh compo- 
nent (centered at u>l), elastically scattered on the dressed 
states I + (N - 1)) and | - (N - 1)). The correspond- 
ing amplitudes add coherently, so that the two scatter- 
ing processes are equivalent to elastic scattering by the 
atomic ground state |1). Since the ground state does not 
have a linewidth, the resulting frequency distribution is 
the same as for the Rayleigh component of the resonance 
fluorescence spectrum, that is, the width 27 of the ex- 
cited atomic levels. In contrast, the processes (d) and 
(j) involve sideband photons centered at either lol + O 
or ujl — O undergo Raman process on | + [N — 1)) or 
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| — (N — 1)), respectively. This leads to an additional 
broadening of the frequency distribution of photons, as- 
sociated with these transitions, by 37. Neither of the 
processes (d) and (j) interferes with (e) and (g), because 
the frequency changes associated with these transitions 
are different, and information about the latter is carried 
into the environment by the undetected photons [l6[ . As 
a result, each of the processes (d) and (j) contributes to 
the CBS ladder spectrum with a Lorentzian centered at 
i^ = 0, with linewidth 67, which adds to the one due to 
the Rayleigh transitions. 

As for the sidebands at v — ±ft, there are three pro- 
cesses contributing to each of them. Since the spectrum 
is symmetric, let us analyse the low-frequency sideband. 
The relevant processes are (a), (c), and (f) in Fig. [7j 
Processes (a) and (c) are completely analogous to (e) 
and (g). That is, they add coherently, and the result- 
ing linewidth is the same as for sideband photons of the 
single-atom spectrum, that is, 37. The process (f) is 
analogous to the process (j), hence the linewidth of the 
respective Lorentzian must be the sum of the linewidths 
of the Rayleigh photon and of a dressed state. Thus, we 
obtain a Lorentzian with linewidth 57 centered at u>l— ft. 

We proceed in the same manner with the sidebands 
at v = ±2fl. The respective diagrams for the lower- 
frequency sideband are (b) and (1). Since, in both cases, 
the sideband photons of the resonance fluorescence spec- 
trum undergo Raman processes, the linewidth of the as- 
sociated Lorentzian is 67. 

Also note that the crossed term has the same lineshape 
as the ladder term, at v = ±2fl [see Fig. HJa)]. This pe- 
culiarity of the CBS spectrum will be discussed in more 
detail in the next subsection. In brief, it is related to 
the fact that there is only one direct and its reciprocal 
process, which contribute to the CBS spectral lines at 
either one of these frequencies. Therefore, these inelastic 
photons (self-)interfere perfectly well. Although the con- 
tribution of the processes at v = ±2fl is very small, it 
is non-negligible. The numbers for the integral contribu- 
tions of the ladder and crossed terms to the resonances 
in Fig. Elja) (these numbers can be easily extracted from 
Eqs. (|B19|) and (|B20j) ) yield precisely the asymptotic 
value of the enhancement factor 0^ = 23/21 derived 
in [23]. Without the contributions from the outer side- 
bands, one would underestimate by approximately 
6%. 

Finally, let us consider the doublet at ojl ± ft/2. The 
respective diagrams in Fig. [Jj (m) and (n), describe 
the resonant scattering of photons on the transitions 
| + (iV-l)) -> |2) and |- (TV- 1)) -> |2). Since the incom- 
ing photons do not originate from the (Mollow) peaks of 
the resonance fluorescence spectrum emitted by the first 
atom, the corresponding linewidths are only defined by 
the linewidths of the 'filtering' transitions of the second 
atom. These linewidths 37 are a sum of the width 27 of 
the excited level |2), and the linewidth 7 of either of the 
dressed states. 



b. Detuned case. For 5^0, the major part of the 
above analysis for 6 — is still valid. Though some 
modifications are needed, in order to explain why some 
resonances creep towards the other ones; their weights 
are redistributed, and even the interference character of 
some of them is changed (see Fig. E£b)). Here, we will 
comment on their linewidths and positions. 

To derive the new linewidths of the ladder term, one 
needs to account for the dependence of linewidths of the 
resonance fluorescence spectrum of the first atom on the 
detuning [2||. With that and Fig. [7J one easily obtains 
the linewidths of the CBS ladder spectrum, at finite S. 

It is easy to show that the new resonance frequencies 
of the CBS transition |1) <-» |2), modified by the ac- 
Stark effect, are w ± (0' ± 5)/2,where ft' = ^ft 2 + 5 2 
is the modified Rabi frequency (26[. It follows that for 
the detuned case, seven CBS resonances manifest at v = 
±2ft' v = ±fl', v = ±(fl' =F S)/2, and v = 0. Therefore, 
as we increase the detuning, the Autler-Townes doublet 
centered at v = —(ft' + 5)/2 approaches the sideband 
at v = —ft', whereas its counterpart centered at v = 
(ft' — 8)/ 2 approaches the central resonance, as evident 
from Fig. Mb)- 

Since the single-atom resonance fluorescence spectrum 
is always symmetric, this also explains why the sideband 
at v — — ft' which is red-shifted with respect to v = 
gains more weight than its blue-shifted counterpart: It is 
closer to the transition frequency ui' a = ui n — {ft' — S)/2 
as compared to the distance of the sideband at v = ft' 
from uj'q = oj a + (ft' — 5)/2. Consequently, the associated 
scattering cross-section is larger. This also explains a 
larger value of the peak at v = —(ft' + S)/2 than at 
v = (ft' — <5)/2, given that, under the detuned driving, 
these resonances originate mainly from the two sidebands 
of the Mollow triplet. However, when we consider Raman 
scattering at v = ±2fl', the scattering cross-section for 
a photon atw^ + O' depends on its detuning (ft' + S)/2 
from the transition |2, N) <-> | + (N - 1)) [see Fig. [7J1)]. 
Likewise, the cross-section for a photon at ojl — ft' is 
determined by its detuning (2>fl' —5) /2 from the transition 
\2,N) <-> I — (JV — 1)) [see diagram (b) on the right of 
Fig. U\. For ft = IOO7 and 6 = 2O7, the latter detuning 
is obviously larger, leading to a redistribution of signal 
weights in favor of the blue-shifted sideband [compare 
the insets of Fig. [5Jb)]. Note also that this asymmetry 
does not spoil the perfectness of the interference of CBS 
photons around v — ±2SY. 



C. Interpretation of the crossed spectrum at 8 = 

Finally, for an intuitive interpretation of the interfer- 
ence character of different components of the crossed 
term's spectrum, we will rely upon the diagrammatic 
technique developed in [jj6[. We will see that a slightly 
modified version thereof in terms of dressed states al- 
lows for the explanation of the destructive interference 
at v = ±ft, of the dispersive resonances at v = ±fl/2, 
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FIG. 9: Qualitative interpretation of the interference char- 
acter. Grey-shaded and black spots represent atoms scatter- 
ing inelastically and quasi-elastically, respectively. Solid and 
dashed arrows correspond to direct and reversed scattering 
amplitudes, respectively. In examples (a) and (b), the direct 
and the reversed amplitudes are reciprocal, whereas they are 
not in cases (c) and (d). 



and of the constructive interferences at the other reso- 
nances. 

Let us first note that the interference character of dif- 
ferent resonances of the crossed spectrum is defined by 
how much, if at all, the outgoing light frequency is shifted 
with respect to the incoming laser frequency ujl- In the 
secular limit 17 ;§> 7, these shifts are of the order of 17, 
with the widths of the shifts' distribution of the order 
of 7. Frequency shifts give rise to phase shifts between 
the interfering amplitudes. Taking into account the con- 
tinuous w-dependence of the amplitudes' phases, we will 
ignore the frequency changes of the order of 7, and refer 
to the respective scattering processes as 'quasi-elastic'. 
We will label 'inelastic' the processes upon which the 
frequency changes by the quantity of the order of 17. 

With these preliminaries, we will use diagrams similar 
to that of Fig. [3J to depict direct and reversed scattering 
amplitudes for the scattering processes contributing to 
the CBS spectrum, in the limit of well-separated spec- 
tral lines (see Fig. [8]). These diagrams, shown in Fig. [9l 
split into three categories: (a), both atoms scatter quasi- 
elastically; (b), (d), both atoms scatter inelastically; (c), 
one of the two atoms scatters inelastically. When both 
atoms scatter quasi-elastically [Fig. [9fa)] , or inelastically 
- into the outer sidebands [Fig. EJb)], there is a recipro- 
cal, reversed process, such that the interference is perfect, 
as indeed observed for peaks at v ~ and v ~ ±217 of 
Fig. Ufa), since there, the crossed and the ladder contri- 
butions arc indistinguishable. However, note that, since 
each inelastic event induces the same change of frequency 
for both reversed processes, the intermediate photons 



(wj,±f2) have opposite detuning from the laser frequency 
in Fig. [H^c). In the case 8 = 0, those two frequencies are 
equally far detuned from the atomic resonance, and, as 
evident from Fig.[5]Ja), the equality between the reversed 
amplitudes remains preserved. Only for non-vanishing 
detuning, the interference contrast at v = is slightly 
reduced, see Fig. [SJb), where the weights of the crossed 
and ladder contribution (0.144 and 0.157) are different. 

Now let us address diagram (d) of Fig. [9] which is rem- 
iniscent of Fig. |U This diagram shows non-reciprocal in- 
terfering amplitudes, which, depending on the frequency 
shifts between the intermediate photons for direct and 
reversed paths, can describe either the dispersive fea- 
tures (which corresponds to v = ±17/2) of the crossed 
term's spectrum, or the destructive interference thereof 
(at v = ±17). In order to see this, we will construct 
explicit expressions for the direct and reversed scatter- 
ing amplitudes. Although both atoms are driven by a 
powerful laser field, it is sufficient to consider the two- 
photon scattering amplitudes, as in Fig. |4j provided that 
the ac-Stark shifts of the CBS transition |2) «-> |1) as 
well as of the laser-driven transition |4) <-» |1) are prop- 
erly accounted for. In other words, for the CBS tran- 
sition, instead of considering the atomic resonance fre- 
quency uj , one should consider two resonance frequen- 
cies oj' = u>o — 17/2 and uj'q = ujq + 17/2. On the other 
hand, ujq remains the resonance frequency for the laser- 
driven transition. It should be stressed that our present 
treatment is very schematic and does not pretend to be a 
full explanation of the interference character in the sat- 
uration regime. Rather, we aim here at a qualitative 
understanding of the interference effect as observed in 
Pig. EI 

In order to construct the direct and reversed ampli- 
tudes, E\ and E2, respectively, we employ the symmetry 
of the CBS spectrum at 5 = with respect to a change of 
the sign of v, and consider only the interference charac- 
ter of the two resonances, at v* ~ 17, 17/2. In the direct 
two-photon process, photons with frequencies ujl + v* 
and u>l — v* are scattered by the first atom. We assume 
that the detected photon with frequency w* = Wl + v* 
is quasi-elastically scattered on an ac-Stark shifted tran- 
sition with frequency uj'q of the second atom, since 
is closer to the modified transition frequency uj'q than to 
the one with frequency uj'q. For the reversed amplitude, a 
laser photon is scattered quasi-elastically by the second 
atom. Then, the first atom scatters inelastically a pair of 
photons with frequency ujl into photons with frequencies 
W* and 2ujl — w*, with the detected photon originating 
from the Raman-anti-Stokes scattering on the ac-Stark 
shifted transition with the resonance frequency uj'q, while 
the undetected photon scatters on the laser-driven tran- 
sition, for which ujq is the resonance frequency. 

With the above assumptions, the amplitudes E\ and 
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E 2 read 



Acknowledgments 



Ei 
E 2 



1 



1 



i<j>/2 



*v* — i"f + Yf) — fi/2 — 17' 
1 1 \ ie**/ 2 



. z/* + fi/2 — 17 is*+ijJ 7 



(23) 
(24) 



It is easy to check that the phase shift between the non- 
reciprocal amplitudes (|23|) and (|24|) continuously depends 
on (or 1/*). For ^„ ~ 0/2, in full analogy with the 
case studied in Sec. IIIIB2I [since v* corresponds to the 
(modified) transition frequency], the interfering ampli- 
tudes are phase shifted by approximately 7r/2, and thus 
describe a dispersive line shape, under slight variation of 
v*. The same amplitudes, for ~ O, are phase shifted 
by approximately 7r and therefore interfere destructively. 
This is quite analogous to the interference-induced anti- 
enhancement of CBS in the regime of elastic scattering, 
due to the non-reciprocity of combined Rayleigh and Ra- 
man scattering processes (33[. Both the dispersive line 
shape and the anti-enhancement is precisely what is ob- 
served in Fig. [51(a). 



IV. CONCLUSION 

We presented a detailed calculation of the spectrum 
of coherent backscattering of light by two identical, ran- 
domly placed atoms, for arbitrary strengths of the laser 
field. 

We saw that the CBS spectrum exhibits its most com- 
plicated structure in the limit of strong driving. The 
ladder term's spectrum reveals seven Lorentzian peaks 
which represent: (i) Mollow triplet photons (re-)scattered 
by the second atom, with resonances at oj = u>i — fi, 
lu = lul, uj = lul + Q] (ii) an Autler-Townes doublet with 
resonances at li = lol — 0/2 and u> = ujl + 0/2; (iii) a 
doublet originating from the Raman Stokes scattering of 
the low-frequency sideband, and from the Raman anti- 
Stokes scattering of the high-frequency sideband of the 
Mollow triplet, with resonance frequencies u> = ljl — 2fi 
and uj = lul + 20, respectively. 

The crossed term's spectrum reveals resonances lo- 
cated at the same frequencies as those of the ladder 
term. We interpreted the interference character of dif- 
ferent resonances through reciprocity arguments applied 
to dressed states. Interference is always constructive for 
the central resonance at u> = ujl, and for the outer side- 
bands at u> = ±217. It changes character in the vicin- 
ity of the Autler-Townes doublet (at u — ujl ± 0/2, 
where the crossed term's spectrum is described by a dis- 
persive curve), and contributes a negative Lorentzian at 

LU = UJL ± f2. 

For the case of the detuned driving, at \8\ — 0, one of 
the dispersive resonances overlaps with one of the neg- 
ative Lorentzians. Interference between the scattering 
processes responsible for appearance of these two reso- 
nances may lead to CBS anti-enhancement. 
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APPENDIX A: CBS INTENSITY AND 
ENHANCEMENT FACTOR 

1. Derivation procedure 

The operator master equation (JSJ) leads to the linear 
matrix equation 



(Q) = (A + V)(Q)+j. 



(Al) 



Here, vector Q having 255 elements is obtained from the 
tensor product qi <8> q2 , where 

q = (l/2,^i/2, jU 2 /2, /x 3 /2,<7i4,o-4i,CTi3,<T3i,eri2, 

021) (T34, 043, O42, (724, C32, 023) T (A2) 

is a vector whose elements comprise the complete or- 
thonormal basis set of operators for a four-level quan- 
tum system, and indices '1' and '2' number atoms. In 

Eq. fijja , 



l = 0"n + C22 

Hi = 022 — 033 
U>2 = 022 — 033 



033 
0"44 
0"44 



/i 3 — 0- 2 2 + 033 — 044 



CT44, 
011, 
011) 

an. 



(A3a) 
(A3b) 
(A3c) 
(A3d) 



In order to describe a mapping of the elements of the two 
16 x 16 vectors onto the vector with 255 elements, it is 
convenient to numerate the basis operators starting from 
'0'. We will denote the ith element of a vector a by [a]». 
For instance, [qi]o = li/2. With these rules, we map 
element ([qi]; ® [q 2 ]m) onto ([Q]«), where n = 161 + m, 
< l,m < 15. After we exclude the first element of the 
tensor product (it corresponds to I = m = 0), we obtain 
that n runs from 1 to 255. Then, the matrices A, V, and 
j are generated by inserting 255 elements of the vector Q 
to Eqs. (O, and performing the quantum mechanical 



averaging: 



25F, 



{(C a + C p )[Q] n ) = A nm ([Q] m ) + [j]„, (A4a) 

255 

{(Cafi + £/Ja)[Q]n) = ^nfn<[Q]tn)- (A4b) 



From the Laplace transform solution to Eq. (|A1|) . 
(Q(z)) = (zt 255 -A- V)- 1 ((Q(0)> + z~% where z 
is defined in the main text after Eq. Q, and I255 is the 
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unit 255 x 255 matrix, we extract the steady state solu- 
tion (Q) ss = lim z ^ z(Q(z)). 

The double scattering contribution corresponds to the 
perturbative expansion of (Q) ss to the second order in the 

coupling constant g, (Q) S s =G VGoVGoj, where G = 
— A -1 . Finally, the correlation functions relevant for the 
evaluation of the total ladder and crossed contributions 
[see Eqs. CCE]), read 



: ([Q]l52)I?, 
([Q]l37)I 2 s ] , 
= <[Q] 

=<[Q]i>£ ] - 



16/ 



H[Q] 

-([Q] 



32, 



.PI. 



(A5a) 
(A5b) 
-([Q]48>I 2 s ] ,(A5c) 



•<[Q]3 



(A5d) 



Concerning the elastic ladder and crossed terms [see 
Eq. (fTS"]) ]. the relevant dipole moment expectation val- 
ues are given by 



[i] 



([q: 



128/ 



[Q]i44)W, (A6a) 



(vW = ([Q] 8 )™, 



■mi 



(A6b) 



where (Q) £ = G VG j. 

Numerical results for the elastic and inelastic intensi- 
ties for different detunings can be found in the main text. 
In the next subsection, we will present analytical results 
for the case of exact resonance. 



2. Analytical results for 5 = 

Evaluating correlation functions from Eq. (|A5|) . we ar- 
rive at following results [20| 



2Re{(a 2 > 12 



<^22>ff + (* 



1 _2 \[2l ik-r 



»} = |g| 2 |A +1 , +1 | 



Ri(s) 



2 \[21 



(4 + s)P(s) 
xcos{(k + k L )-n 2 }, (A7) 



l<?l 2 |A +1 ,+n 



P(s) 



(A8) 



i?i(s), i?2(s)i and P(s) are polynomial expressions in the 
on-resonance saturation parameter s = £l 2 /2'y 2 , 



Ri(s) 
P(s) 



(6912s + 3168s 2 
264s 3 + 20s 4 + s 5 ) 



(A9a) 



^1152s + 528s 2 + 132s 3 + 7s 4 ) , (A9b) 



= (l + s) 2 (12 + s)(32 + 20s + s 2 ), (A9c) 



and A + i i+ i = e +1 • A • e +1 . 

The configuration average of (|A7|) and (|A8|) leads to 
the final result 



cT(0) 



T tot 

L 2 



I + s)P(s)V 15 35 J' (A1U) 

(All) 



(4 + s)P(s) 

■2\g\ 2 R2(s) 



15P(s) 



with g = g\ Ta „ = £. The scattering angle = 2 arcsin{|k + 
ki|/2fcL} "C 1 with respect to the backscattering direc- 
tion was assumed to be sufficiently small herein. 

The enhancement factor a(s), Eq. (TT3")) . deduced from 
Eqs. (|AT0)> and (fATT|) reads 



a(s) = 1 



Ri(s) 



(4 + s)P 2 (s) 



(A12) 



and a(0) = 2.0 in the weak field limit, as it must be in 
the elastic scattering regime. For small s, enhancement 
linearly decreases as 2 — s/4, in full agreement with the 
diagrammatic theoretical result [l6j and in qualitative 
agreement with the result of Sr experiment 8]. When 
s increases further, a monotonically drops to an asymp- 
totic value lim^oo a(s) — — 23/21 [20j | which is 
strictly larger than unity, implying a residual construc- 
tive interference in the deep saturation regime. 

We will next show that this interference is due to in- 
elastic photons only. Indeed, we obtained the following 
result for the elastic ladder and crossed terms 



p; 



cf 



15 (1 + s) 4 



(A13) 



As seen from Eq. (|A13|) the elastic component shows per- 
fect contrast for all s. In particular, it is this component 
that results in enhancement a = 2 for very small s — > 0. 
However, in the deep saturation regime, this component 
decreases as s~ 3 , while the counterparts of the total in- 
tensity, Eqs. (|A10|) . (|A11|) . as s _1 . Herefrom follows our 
conclusion about the origin of the residual enhancement 
in the deep saturation regime. Explicitly, the inelastic 
crossed and ladder terms obtained by elementary sub- 
straction of Eq. (|A13[) from Eqs. (|A10j) and (jAllj) read 



£*mcl 



2|g| 2 20736s 2 



2s 7 



L 



incl 



15 9(1 + s] 
2\g\ 2 2016s 2 



»(4- 



s)P(s) 
f 7s 6 



15 



3(1 + s) 2 P(s) 



(A14) 
(A15) 



where we have retained only the lowest and highest order 
terms in the numerators, because these terms will be used 
to verify the expressions for the spectra. Explicit val- 
ues of other coefficients are not important. Using (|A14[) 
and (|A15|) , it is easy to verify that lim s ^ 
2/21 = aoo - 1. 



£fmel J jjiricl 



APPENDIX B: DERIVATION OF CBS 
SPECTRUM FROM EQ. © 

1. Elastic and inelastic spectra 



By virtue of the quantum regression theorem [24| , the 
temporal correlation functions of a Markov process obey 
the same equation of motion as the expectation values 
of operators, that is, Eq. ©, but with different initial 
conditions and a different free term. The latter can be 
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straightforwardly determined provided that the station- 
ary solution to Eq. ([5]) is known. 

The four correlation functions appearing in the defini- 
tion of the spectrum [see Eqs. ([3]), (@|] can be extracted 
from the following correlation functions: 



s a (t) = <o-§iQ(t)}« (a = 1,2), 



(Bl) 



where s a (t) is a vector containing, like (Q(t)), 255 ele- 
ments. The vector s a (t) satisfies the equation of motion 



S a = (A + V)S Q + (cr 21 ) ss j, 



(B2) 



which explicitly accounts for the modification of the free 
term. From the definition (IBlj) , it follows that the vectors 
si and s 2 are obtained from (q' : (gi q2) ss and (qi ® q 2 ) ss , 
respectively, where 

q' = (o-2l/2,-Cr 2 l/2,0-2l/2,-(72l/2, 0-24,0,^23,0, 

a 22 ,0,0,0,0,0,0,0) T , (B3) 

by using the same mapping rule as described previously 
in Sec. lA 11 Finally, to find the initial conditions for these 
vectors, one just needs to make certain transpositions 
among the elements of the vector (Q) ss - 

The Laplace transform solution of Eq. (|B2|) reads: 



s a (z) 



1 



zt 



255 



(A + V) 



MO) + -(^ Q i> 



(B4) 



where s Q (z) is a Laplace image of s Q (t) . Using this solu- 
tion and definitions 



'12/ss 



= 2([qi] 8 <g> [q 2 ] ), (cr 2 2 ) ss = 2([qi] ® [q 2 ]s) 



we arrive at the following expression for the Laplace im- 
age of the correlation function 

GH\z) = 2([Si(*)] 128 + [ia(z)] 8 

+ [s 1 (z)] 8 e 4k - ri2 + [i 2 (z)] 128 e- lk - ri2 ),(B6) 

The next step will be to extract the lowest-order in the 
coupling g, nonvanishing contribution from (|B4|1 and in- 
sert it into Eq. (|B6|) . This contribution is on the order g 2 
and corresponds to double scattering in which the atoms 
exchange a photon. Explicitly, 



s™(z) = Go(z)VG (z)VG (z) (0) + ± (ofr 



Go(z)VG (z) ^(0)+ J -(^ 1 )g 



G (z) 8^(0) + ^ 



J 2 ] 



JO] 



(B7) 



where Gq(z) = (zfl-255 — -A-) -1 , and the superscripts [0], 
[1], and [2] indicate terms on the order g°, g l , and g 2 , 
respectively. 

Expression (|B7[) can be simplified since (|B6[) gives a 
spectrum in the h \\ h channel. Therefore, the station- 
ary expectation values of operators related to the laser- 
nondriven transition vanish in g°. So, the first, second, 



and last terms must be dropped from (|B7|) . Concern- 
ing the latter term, it does not contribute to (|B6|) since 

[Go(^)j]s and [G (z)j]i28 give (cr 2 2 )Is and (ctJ 2 )S, re- 
spectively, and thus vanish for the reason already indi- 
cated. 

It is useful to split the nonvanishing part of (|B7|) into 
two counterparts from which the elastic and inelastic 
spectra are extracted: 



= sL(* > ,<•->• 

where the two vectors 



S [2] 



€LiW = Go(z)VGo(z)sW(0) + Go(z)sL 2 ](0) 



1 



(B8) 
(B9a) 



[G„(z)VGo(z)-G VGo]j, 



(B9b) 



lead to the elastic and inelastic spectra, respectively. 
Inserting (|B9aj) into (jB6]) we obtain 



1 \[i] 

7 21/ss \ <J 12/ss 



-2Re{(a 2 i 1 )W(af 



1 \[l]/^2 ik.r 12 



} , (BIO) 



with the right hand side expression in the round brack- 
ets being nothing but the stationary elastic intensity If 1 , 
Eq. {15|). Finally, putting [Git (z)] e i to (@| we arrive at 



where we have used the formula 
lim 



r^o 1 — iv v 



(Bll) 
(B12) 



with P denoting the principal value of an integral. 

One can check that the right hand side of Eq. (|B9b[) 
does not have a pole at z = 0, which means that the 
respective expression describes the inelastic spectrum. 



2. Analytical results for inelastic spectrum at 8 = 

a. Weak field (£1 < >y) 

At small Rabi frequencies, the analytical formulas are 
obtained after taking into account the lowest-order in- 
elastic process - the two-photon scattering, - and ne- 
glecting the inelastic processes of higher orders. The 
two-photon processes are proportional to the square of 
the intensity, that is, to f2 4 . The ladder and crossed 
terms read (we omit the common prefactor 2|g| 2 /15): 



n\ 4 7 3 (2 7 2 + j,2) 



7 J 2( 7 2 + 1/2)3 



7 J (7 2 + I/ 2 ) 3 



(B13) 
(B14) 
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It is easy to check that the expressions (|B13|) . (|B14I) are 
consistent with the behavior of the enhancement factor 
in the two-photon scattering regime. Integrating {v) , 
C™ cX (v) over all frequencies, we obtain the following in- 
elastic ladder and crossed terms for small il: 



t inel 

L 2 — 



^fincl 



7 fn 



16 V 7 

3 fn 



(B15) 
(B16) 



Combining Eqs. (|B15[) . (|B16[) with the small- s expression 
for the elastic ladder and crossed terms Lf = C^ 1 = s, 
and rewriting Eq. (|B16|1 in terms of s, we recover the 
expected linear decrease 



1 



3s 2 /2 s 
■7s 2 /4~ 4' 



b. Strong field (fl S> ~f} 

In the opposite limit of a strong field, the CBS inten- 
sity is inversely proportional to the laser field intensity. 
We will now present the analytical expressions for the 
ladder and crossed spectra derived in the leading order 
~(7/«) 2 . 

In this case, explicit expressions for CBS spectra can 
be represented by using a function of two real variables 
Xi and X2- 



£{xi,x 2 ) = — 



7r xi + xi 



(B18) 



Let us mention the properties of £{x\,X2) that are im- 
portant to us: (i) if X\ — Const, then the function (|B18p 
represents a Lorentzian with full width at half maximum 
(referred to as width in the main text) 2x\ and resonance 



at X2 — 0; (ii) if x 2 — Const, then (|B18[) describes a 
resonance of a dispersive type at x\ = 0, with the width 
2cc 2 . 

With the help of the function (|B18[) . the ladder and 
crossed spectra are given by 



~LT\v) 
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+y [^(37/2, v - fi/2) + £(3 7 /2, v + n/2)] 
+~[£(3 7 /2, v - 0) + £fry/2, v + n)] 
+ ^[£(57/2, u-Sl) + £(5j/2, v + fi)] 



-—[£(*y,v-2tl) + £(*y,v 



20)] ,(B19) 



(B17) CT» 



6 



-£{2 1 ,u) + -£{^,u) 



[£(5 7 /2, v - O) + £(57/2, v + ii)} 



+ ^[£(3 7 , v - 2fl) + £(3 7 , v + 2fi)]) 
/7\3 208 

+ \n) 45" 

x [£{v + n/2, 37/2) - £(1/ - 0/2, 3 7 /2)], 

(B20) 

where the two terms of order (7/T2) 3 are retained be- 
cause they define dispersive resonances of Ci, 1101 ^) at 
v = ±0/2. By performing the elementary integrations 
of Eqs. (|B19|) and (|B20|) we arrive at the inelastic ladder 
and crossed terms 



j mcl 
L 2 



14 (1 
~3 



Ci 1 



(B21) 



which are consistent with Eqs. (|A14[) , (|A15|) and, hence, 



with a = otr 



23/21. 
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